Spin relaxation in narrow wires of a two-dimensional electron gas 
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How does an initially homogeneous spin-polarization in a confined two-dimensional electron gas 
with Rashba spin-orbit coupling evolve in time? How does the relaxation time depend on system 
size? We study these questions for systems of a size that is much larger than the Fermi wavelength, 
but comparable and even shorter than the spin relaxation length. Depending on the confinement 
spin-relaxation may become faster or slower than in the bulk. An initially homogeneously polarized 
spin system evolves into a spiral pattern. 
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INTRODUCTION 



Long spin relaxation times are desirable for the opera- 
tion of spintronic devices. In confined systems like quan- 
tum dots with a discrete level spectrum an initial electron 
spin polarization will not decay exponentially, unless the 
spins are coupled to extra degrees of freedom like phonons 
or nuclear spinsi^ For larger or open systems, where the 
energy spectrum can be considered continuous, exponen- 
tial decay is possible by purely elastic scattering which is 
then the dominant process at low enough temperatures. 

In this article we reconsider spin-relaxation in a two- 
dimensional electron gas confined to a narrow channel, 
concentrating on elastic scattering processes. Our study 
has been motivated by a recent experiment, where the 
size dependence of the spin relaxation time of an n- 
InGaAs wire has been measured via Faraday rotation 
spectroscopy^ The size dependence sets in when the 
wire width is of the order of several bulk spin relax- 
ation lengths. The observed behavior is non-monotonic 
with an initial increase followed by a sharp decrease at 
the smallest wire widths. The maximum spin relaxation 
time occurs when the width is of the order of the bulk 
spin relaxation length. Remarkably, while the tendency 
to suppress spin-relaxation in confined systems has been 
predicted in a number of theoretical works there 
has been no anticipation of the increase of the spin re- 
laxation observed at the smallest widths. In this article 
we show that spin active boundaries, not considered in 
the previous theoretical analysis, dramatically change the 
size dependence of the spin relaxation time in the small 
width limit and provide a useful point of view as far as 
the interpretation of the experiment is concerned. 



The rest of the paper is organised as follows. In Sec. 
II we recall the method of quasiclassical Green functions 
and point out when a diffusion equation approach is pos- 
sible. In Sec. Ill we consider the spin dynamics in the 
bulk, whereas in Sec. IV we turn to spin relaxation in 
finite systems. A summary is given in Sec. V. 



II. THE BASIC EQUATIONS 

We start from the Hamiltonian of a two-dimensional 
electron gas 



H = i_+b.<T + V(x), 
2m 



(1) 



where b is the internal magnetic field due to the spin- 
orbit coupling, er is the vector of Pauli matrices and 
V(x) is a delta-correlated random potential due to im- 
purity scattering. For the internal field we concentrate 
on the Rashba model, b = ap x e 2 , where the spin- 
orbit field arises due to a structural inversion asymme- 
try. However, our treatment can be applied also to 
the case of bulk inversion asymmetry, where there is 
in addition the Dresselhaus termifi which contributes 
to the two-dimensional effective Hamiltonian with b = 
P(p x ,—Py) + j{PxPl,-p y pi), or to the case of a two- 
dimensional hole gas, where the Rashba spin-orbit field 
has a cubic dependence on momentum^* 

To study the spin dynamics we rely on the method of 
the quasiclassical Green function»i2ii^ The latter solves 
the Eilenberger equation (h= 1), 



d t g + v F - Vg + i[b ■ cr,g] 
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(g -(g)), (2) 



where r is the elastic scattering time, arising by the 
adoption of the standard self-consistent Born approxi- 
mation for elastic impurity scattering. v> is the Fermi 
velocity and (...) denotes the angular average over the 
Fermi surface. Eq. © is valid in the limit when both 
the spin-orbit energy and t -1 are small compared to 
the Fermi energy^* Besides the space and time depen- 
dence implied by Eq. (J2J, the quasiclassical Green func- 
tion g = <7 ss '(e, Vpi x, t) is a matrix in spin space and 
is a function of energy e (measured with respect to the 
Fermi surface) and direction on the Fermi surface, 
In the absence of external potentials and magnetic fields 
the particle and spin density are obtained by summing 



2 



over energy, angle and spin 



p(x, t) 



s(x,t) 



de )(g ss (e, v F ;x,i)), 



(3) 



>(g s > s (e,v F ,x,t)). (4) 



When the spin-orbit energy is small compared to the 
scattering rate the spin-dynamics becomes diffusive in 
the low frequency, long wavelength limit. In this limit, 
all the harmonics characterizing the angle dependence of 
g can be expressed conveniently in terms of the s-wave 
component. As a result the problem simplifies consider- 
ably and one obtains, for the angle averaged Green func- 
tion, (g), a diffusion equation, which, in the presence of 
a general spin-orbit field, acquires the form given, for 
instance, in Ref. For the Rashba Hamiltonian in par- 
ticular the equations read 
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2Cd x s z 



2Cd y s z 



2Cd x 



(5) 
(6) 

(7) 



2CdySy, (8) 



Since here we focus on the spin dynamics, we limit our- 
selves to terms to leading order in the parameter a/vp, 
therefore neglecting spin-charge coupling 1416 ' 17,18 . 



III. SPIN DYNAMICS IN THE BULK 

The spin-dynamics is particularly simple for a spatially 
homogeneous spin density. From the diffusion equation 
one observes that the different spin directions decouple 
and the spin polarization along the e XtV or e z axis re- 
laxes with the time constant t s or r s /2, respectively. This 
type of spin relaxation, known as the D'yakonov-Perel' 
mechanismjSi is easily understood: the electron spin pro- 
cesses around the internal field b. Scattering from an 
impurity changes the direction of the internal field, and 
thus randomizes the spin precession. If the precession 
angle, (5$ ~ |b|r, between collisions is small, the spin 
dynamics resembles a diffusion. The total precession an- 
gle after TV scattering events is of the order \/N53>, from 
which the spin relaxation rate is estimated of the order 
t s ~ t/8$ 2 ~ T/(ap F T) 2 . 

In general we cannot rely on the diffusion equation but 
have to solve Eq. In the spatially homogeneous case 
the task simplifies, since again different spin directions do 
not couple. As a result, the spin dynamics is described 



L(L 
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L 2 - -a 2 
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[L(L 2 



L 2 s, 



(9) 
(10) 



with L = 1 + rdt and a = 2api?r. Apparently the 
full description of the spin dynamics now requires three 
time constants for each component. The slowest decaying 
component goes as Si(t) ~ exp(— jit) with 



Jx,y 
Jx,y 

Iz 



1/Ts (a < 1) 
l/2r (a > 1) 



(11) 

(12) 
(13) 



i.e. for a clean system all components decay on the time 
scale of the scattering time r. 

For conventional diffusion the mode with the longest 
lifetime is homogeneous in space. Here we encounter a 
different situation, due to the coupling of the various spin 
components. From the diffusion equation, for instance, 
the modes with the longest lifetime are found to form an 
elliptically modulated spin spiral, 

s(x, t) oc [e q cos(q • x) + Ae z sin(q • x)] exp(— jot), (14) 

where the vector q lives in the x-y-plane with Dq 2 = 
15/16t s , A — 3/y/TE, and 70 = 7/16r s . These long-living 
modes have been investigated in and in p5T |. 



IV. SPIN DYNAMICS IN FINITE SYSTEMS 

For finite systems the equations have to be supple- 
mented by boundary conditions. 22 ' 23 24 For specular scat- 
tering where an in-going trajectory is scattered into one 
outgoing direction the boundary condition for the Eilen- 
berger equation reads^ 



fJ, 



(15) 



The unitary matrix S describes the surface scattering. 
By decomposing the Green function in charge and spin 
components, g ss > = go5 SS ' + g ■ Css' , the boundary condi- 
tion is equivalently expressed as 



out 

9o 



96 



(16) 



with an orthogonal matrix R that rotates the spin at 
the surface. Charge conservation implies that no current 
flows through the boundary, 



(n- VFffo) oc n-j c = 0, 



(17) 



where n is a vector normal to the boundary. For a spin 
conserving boundary (R — 1) also all components a of 
the spin current perpendicular to the surface are zero, 



(n • VFg a ) oc n ■ j a = 0. 



(18) 



For the general case, R ^ 1, Eq. (|18|) is not valid. 

In the following we consider two types of boundary 
conditions. First we consider a spin-conserving bound- 
ary, where 



Ikin S) 



|kout s), 



(19) 
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such that S and R are unit matrices. As a second exam- 
ple we will consider a boundary that scatters adiabati- 
cally, i.e., an incoming wave in an eigenstate |ki n ±) of the 
Hamiltonian - not including disorder - is scattered 
into the same band, 



|k in ±) 



(20) 



as it is expected for a smooth confining potentialiS*^ 
Since the eigenstates in the presence of the field b are 
spin-polarized parallel to b, the adiabatic boundary im- 
plies a rotation of the spin polarization, S = e 10zV , where 
(p is the angle between the electron momentum and the 
normal of the interface. 

Before presenting our numerical results for the spin re- 
laxation obtained by solving Eq. Q , it is useful to derive 
the boundary conditions in the diffusive regime. The idea 
is to solve the Eilenberger equation near the boundary as- 
suming that the angular average of g varies only slowly 
on the scale of the mean free path, I. This is justified 
when the spin relaxation length L s is much longer than 
I. For directions pointing into the boundary one can then 
use the expansion valid in the bulk, 



(g) - rv F ■ V (g) + 2rb x (g) 



(21) 



Combining i|16|l and i|21|) yields a relation that connects 
linearly the three spin components of the angular aver- 
aged Green function, (g), and its spatial derivatives. The 
resulting form of the boundary condition depends on the 
rotation matrix R(<p), 



2^ 



[1 + R] [(g) - rv F ■ V (g) + 2rb x (g)] = (g) 



(22) 

For the Rashba model, using the expression for the charge 
and the spin current in the diffusive limit, we get for spin- 
conserving boundary conditions'^ 



Dd x s x - Cs z = n • jj, = 0, 
-Dd x s y = n ■ jy = 0, 
-Dd x s z + Cs x = n j 2 = 0, 



(23) 
(24) 
(25) 



where n is in the ^-direction. For adiabatic boundary 
conditions, in contrast, angular averaging of Eq. 122|) 
yields that 



and s y = 0, 



(26) 



while the z-component of the spin is still conserved and 
therefore Eq. (|25|l remains valid. 

Fig. ^ shows the time evolution of the spin profile for 
a long wire of width L = 2001, where I = vft is the elas- 
tic mean free path. Initially the spin was homogeneously 
polarized in z-dircction. The results were obtained from 
the Eilenberger equation with ap F r = 0.1, and the con- 
serving boundary condition, S ss i = 6 ss i. Inside the wire 
one observes a homogeneous decay of the spin polariza- 
tion, with the time constant r s /2. At the boundaries long 
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FIG. 1: Time evolution of the spin polarization in a wide 
channel (L — 200 1 m 40 L s , dp ft = 0.1). The curves from 
top to bottom correspond to different times, with At = 50r sa 
t s . s z changes sign at various positions where a steep drop of 
\s z \ is visible in the figure. 




FIG. 2: Lowest eigenvalues of the spin-diffusion operator 
for the Rashba model and conserving boundary conditions, 
Eqs. 1231 — (12511 . Modes with 7 < 7/16r s have a complex wave 
vector, and can therefore exist only at the edges of the wire. 
The dashed curve is 7t s = (L/L 3 ) 2 /12 obtained in Ref.Q for 
very narrow wires. 



living modes show up which dominate the spin profile in 
the long time limit. 

For further investigation of these modes we write the 
spin diffusion equation, Eqs. © ©■ as dts + js = 0, and 
determine the eigenvalues and eigenmodes of the opera- 
tor 7. The eigenmodes are superpositions of plane waves. 
The low frequency spectrum of 7 is shown in Fig. [21 as 
a function of the wire width. Recall that the smallest 
eigenvalue for a bulk system is 70 = 7/16t s . The modes 
with smaller decay rate have a complex wave vector and 
are thus localized near the edges of the wire. For a wide 
system we find a continuum of eigenvalues above 70, and 
two localized modes at 7 w 0.382/r s . 

For a narrow wire most strikingly one eigenvalue goes 
to zero with decreasing width, asymptotically as 



1T S 



12 V L s 



(27) 



This corresponds to the suppression of spin-relaxation in 
small systems reported earlier by other authorsA^LS^ 
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This effect can be traced back to the specific form of 
the spin-orbit field in the Rashba Hamiltonian being 
proportional to the velocity^ Here we formulate the ar- 
gument for a system including Rashba and also linear 
Dresselhaus term within the spin-diffusion equation ap- 
proach. For a spin profile that is homogeneous in the 
y-direction the angular averaged Eilenberger equation, 
Eq. 0, yields 

d t s x + d x f x = -2maf z -2mPj y z (28) 
d t s y + d x j* = -2maj z >-2m(3f z (29) 
dt-s z + <hji = 2ma(f x + f y ) + 2mp{j v x + f y ) (30) 

where j£ is the spin current polarized in a and flowing 
in /3-direction. In the diffusive limit the spin current 
densities are given by 



J a 



-Ddf3S a + 2r{v P F {h x s) Q ) 



(31) 



which allows to reproduce the spin diffusion equation, 
Eqs. 111 narrow systems the slow modes have 

a smooth density profile, such that to leading order in 
the system size the current can be considered constant 
in space. For a quantum dot, i.e. a system that is con- 
fined in all spatial directions, the vanishing of the spin 
current through the boundaries then implies immediately 
that dts = 0. For a narrow wire the situation becomes 
slightly more complicated since only currents flowing into 
the boundary are zero, which after some algebra leads to 



d t s 



1 



(3 2 a/3 
af3 a 2 





(32) 



In the absence of the Dresselhaus term (/3 = 0) this means 
that d t s x = 0, and d t s v = -s v /t s , d t s z = —s z /t s , which 
implies that the long-living mode in Fig. [3 is polarized 
in x-direction. In the presence of both a Rashba and a 
Dresselhaus term, the spin is still conserved for one di- 
rection which depends on the relative strength of the two 
terms: Perpendicular to the boundary when the Rashba 
term dominates, parallel to the boundary when the Dres- 
selhaus term is larger, and somewhere in between (but al- 
ways in-plane) when both terms are comparable in size. 

The above results change considerably when different 
boundary conditions are applied. Fig. [3] shows the time 
evolution of a spin polarization using adiabatic boundary 
conditions, Eq. I|26|l . Here the spin has been prepared in 
x-direction, i.e. perpendicular to the boundary. In this 
case the boundary mode is absent, and the asymptotic 
decay of the spin polarization is ruled by an inhomogc- 
neous, but extended mode. The spectrum of the spin 
diffusion operator is shown in Fig. ^ The boundary con- 
dition implies that the eigenmodes are s x _ y oc sm(gir), 
s z oc cos(qx) with q — nir/L; the eigenvalues are given 
by 
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FIG. 3: Time evolution of the spin polarization in x-direction 
for a wire with adiabatic boundary condition for the same set 
of parameters as in Fig. Q 




7 (g) = Dq 2 + 



2r, 



2r s V 



16L 2 q 2 



(33) 



FIG. 4: Lowest eigenvalues of the spin-diffusion operator for 
adiabatic boundary conditions 



By inserting the allowed q- values the spectrum shown in 
Fig. 0] is reproduced. In contrast to the previous case of 
a spin-conserving boundary, here all the diffusion modes 
show an increasing spin relaxation rate at the smallest 
wire widths. In particular, all the modes show a non 
monotonous behavior as a function of the wire width with 
a minimum at L/L s — (An/ylo) n m 3 n, where n is the 
mode index. For width L < L s all spin diffusion modes 
relax fast. 



V. SUMMARY 

Motivated by a recent experiment^ we studied spin re- 
laxation in narrow wires with spin-orbit coupling. To 
this end we have solved the Eilenberger equation and 
the spin diffusion equation in the presence of spin-orbit 
interactionia supplemented by boundary conditions for 
both charge and spin degrees of freedom assuming trans- 
lational invariance along the wire. 

The spin diffusion equation has extra, off-diagonal gra- 
dient terms that couple the different spin directions. Due 
to these terms there exist inhomogeneous spin-profiles, 
typically spin-spirals, that decay slower than a homoge- 
neous spin configuration and that show up close to the 
edges of the wire. Therefore there is a tendency for the 
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spin relaxation to slow down, the stronger the influence 
of the boundary becomes, i.e. when decreasing the width 
of the wire. 

For very narrow wires, i.e. narrower than the spin re- 
laxation length, the spin relaxation rate depends crucially 
on the boundary condition and on the form of the spin- 
orbit coupling. For models with linear-in-momentum 
spin orbit field and spin-conserving boundaries the spin 
relaxation rate goes to zero with the system size. It is 
interesting to note that while in a quantum dot all spin- 
components are conserved in the finite size limit, in the 



case of a narrow wire only one of them is conserved. The 
direction of this conserved component depends on the 
relative strength of the Rashba and Dresselhaus terms. 

On the other hand in the case of spin-active bound- 
aries, the spin is not conserved and the relaxation rate 
grows when the wire becomes very narrow. We believe 
that the experimentally observed^ sharp increase of the 
spin relaxation rate in very narrow wires can be explained 
by spin scattering at the boundary. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft through SFB484. 
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